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Photon bunching in parametric down-conversion with continuous wave excitation 
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The first direct measurement of photon bunching (<j( 2 ) correlation function) in one output arm 
of a spontaneous-parametric-down-conversion source operated with a continuous pump laser in the 
single-photon regime is demonstrated. The result is in agreement with the statistics of a ther- 
mal field of the same coherence length, and shows the feasibility of investigating photon statistics 
with compact cw-pumped sources. Implications for entanglement-based quantum cryptography are 
discussed. 
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I. INTRODUCTION 

Light sources based on spontaneous parametric down- 
conversion (SPDC) in which photons from a pump 
laser are sporadically converted into pairs of photons at 
lower frequency in a non-linear crystal, have become an 
essential component in the toolbox of quantum optics 
laboratories in the last decade. These sources allow state- 
of-the-art experiments on entanglement and quantum in- 
formation, due to the ability to supply photons that are 
highly correlated in time, energy and polarisation. 

SPDC-based systems have been exploited as a basis for 
entangled-photons schemes (H as well as approximations 
of single-photon sources in the so-called heralded-photon 
schemes |3(. Such setups are however limited by the 
possibility that additional pairs are created within the 
window determined by the pump-pulse duration or by 
the electronic gate width of the single photon detectors 
(for continuous sources). In multi-photon interference 
experiments, additional pairs in general affect the purity 
of the investigated state In quantum cryptography 
H instead, multi-photon pairs may constitute a security 
threat, since their information-carrying degree of freedom 
(typically polarisation) might be correlated 0]. 

The topic of SPDC-pair emission and correlation, 
deeply related to the statistical properties of the down- 
converted field, is therefore not only interesting for fun- 
damental quantum optics experiments, but also concerns 
applied quantum information. Significant theoretical and 
experimental investigations have already been conducted 
on this subject, but a direct demonstration that the field 
of one arm of an SPDC source is thermal is still missing 
for the special case of a source in single-photon regime 
operated with a continuous pump. Such evidence is pre- 
sented in this communication by measuring the second- 



order correlation function with a Hanbury Brown- Twiss 
setup 0] and comparing it to a prediction based both on 
the detector jitter and a first-order correlation measure- 
ment with a Michelson interferometer. This approach re- 
lies entirely on experimentally observable quantities and 
no assumptions are needed. 

The intrinsic difficulty of the experiment stems from 
the fact that resolving times are much larger than the 
coherence time t c of the down-converted field, so that 
all relevant quantities are averaged out and the statis- 
tics becomes similar to a Poissonian one, with the true 
(thermal) statistics often overlooked. 

The remaining part of the introduction presents a re- 
view of the theoretical background (Sec ll A[) and exist- 
ing experimental evidence (Sec ll B[) concerning multi-pair 
emission, as well as a discussion of the potential security 
risk for quantum cryptography (Sec ll C[) . The experi- 
mental setup, data acquisition (Scc|TT|, and quantitative 
expectation model (Sec lIII[) are then described, followed 
by a discussion of experimental findings (Scc lIV [) . 



A. Field statistics and theoretical description of 
SPDC 

A direct characterisation of the temporal statistical 
properties of a generic quantum electromagnetic field is 
provided by the (normalised) correlation functions, de- 
fined as [1, ch.6] 

(H ' (EHt)E(t)) G (1) (t,t) { ' 



5(2) (M 



(£t(t) E\t + r) E(t + t) E(t)) 
(&(t)E(t)) (ft(t + T) E(t + r)) 

G (2) (t,t + T) 
G {1) (t,t)G {1) (t+T,t+TY 



(lb) 

(1c) 



* These authors contributed equally to the work, 
t Corresponding author: hannes.hucbcl@gmail.com 



where E are field amplitude operators. These quantities 
are, respectively, field and intensity correlations at times 
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t and t + t at the same spatial location. Since for large 
time delays realistic fields are uncorrelated, 3(2) (oo) = 1. 
The behaviour for finite delays depends however on the 
actual statistics of the field. If <?( 2 ) increases around r = 
0, the field is said to be bunched. 

A chaotic field in an interval short with respect to its 
coherence time r c will show thermal statistics, that is the 
distribution P n of the number of photons is 

v n 

Pn = + (2) 

where v is the average number of photons in the interval. 
Chaotic light, for which gr 2 ) (0) = 2, is therefore bunched. 
In contrast, the output of an ideal mono-mode laser dis- 
plays Poissonian behaviour, that is g(2)( T ) — 1- Hence, 
in a short time interval, the probability of counting two 
photons for chaotic light is twice as large as that expected 
for Poissonian statistics. 

With the picture of an SPDC process as a spontaneous 
emission of pairs from the splitting of pump photons, it 
could be naively believed that the statistics of pairs is 
similar to that of the pump field. It is however known 
Q that linear interaction in a resonant medium affects 
incident light in such a way that it becomes a mixture of 
the initial "signal" field and an additional thermal field 
from spontaneous emission. If the incident field is absent 
(or already thermal) , the output field is therefore exactly 
thermal. In analogy, since the SPDC process corresponds 
to optical amplification without input signal (as both sig- 
nal and idler fields are initially in their vacuum states), 
the statistics of pairs needs not mimic that of the pump, 
and thermal behaviour is to be expected. 

As for incoherent light sources, the quantum inter- 
pretation of SPDC bunching relies on constructive two- 
photon interference between probability am plit udes for 
the paths of two indistinguishable photons [ToL [ill ]. It 
has been shown that this process can equivalently be in- 
terpreted as a pair emission stimulated by another down- 
converted pair [ll], ch.7]. 

That the statistics of one arm of a mono-mode SPDC 
source is thermal was theoretically first proven in 1987 
by Yurkc and Potasek [T2|, who derived from the pure 
state of the combined single-mode down-converted fields 
the density operator of the mixed state of one field only 
by tracing out either signal or idler field. In this model, 
the quantum state \<j>) of the combined signal-idler fields 
in the Fock state expansion reads 

oo / \ n 

e i (j)a t 6 t +7) . o6) = y> I r, \ 

^ VH J 

= (f-ib ? | 2 )|0)+r ? |l) si +^|2) sl + O(H 3 ), (3) 

where 77 is a parameter proportional to the pump ampli- 
tude, a) (6 t ) is the creation operator for the idler (signal) 
mode with fixed polarisation, and \n) si = -=j (er&i)"|0) is 
a Fock state with n pairs. By tracing over one field, say i, 
the dominant signal-idler correlation is suppressed, and 



a mixed density matrix is obtained for the residual field 
that displays exactly thermal statistics: 

Ps = J2 p n\n)s(n\ = J2 ^ + 1) (n+l) l")sH, (4) 
n=0 n=0 ^ ' 

where v = sinh 2 177I and P n is the probability of finding 
exactly n signal photons. The SPDC process is neverthe- 
less inherently quantum mechanical, and its "thermality" 
is not an effect of the trace operation, but originates di- 
rectly from the nature of the pair production process [l2| . 

Since SPDC production is far from being mono-mode, 
Yurkc and Potasck's analysis was not realistic; however, 
the thermal nature of a multi-mode SPDC process was 
later confirmed by Tapster and Rarity II 311 , and in a more 
general way by Ou, Rhee and Wang [J, , who analysed 
the case of sources with pulsed pumping. In their deriva- 
tion, what is actually measured is the pulse integral of 
Gm and G( 2 ), i-e., the probabilities of one or two pho- 
tons in the pulse. With time resolution limited to the 
whole pulse, one must resort to test the "bunching ex- 
cess" , which was shown Q to be essentially the ratio of 
the coherence time t c of the photons in the inspected arm 
and the duration t p of the pulse, when t c <C t p . Shorter 
pump pulses then increase the excess, but, contrary to 
intuition, the limit of a very short pulse is not sufficient 
alone to guarantee a full bunching peak (3(2) (0) = 2). 
In (TTL ch.7] it is shown that this condition can only be 
achieved with the help of narrow spectral filtering. 

SPDC sources have also been implemented with con- 
tinuous pumping, using continuous wave (cw) lasers, 
which arc both more practical and more affordable than 
femto-second (fs) pulsed SPDC sources. The field emit- 
ted by cw sources, and its statistical properties, is still 
relatively unexplored. Theoretically, it is usually as- 
sumed that resolving times are so much larger than the 
coherence time of the down-converted field that all rel- 
evant quantities are not accessible. It has been stated 
[H, G3, [H G3 that the statistics of SPDC as a multi- 
mode process, i.e. when the pulse duration of the laser 
is much longer than the coherence time of the produced 
photons, will lead to Poissonian statistics. Careful exper- 
imental design and data analysis, as shown in this paper, 
reveal however the inherent thermal statistics, hidden by 
an apparent Poissonian behaviour. 

B. Current experimental evidence 

Photon statistics in SPDC processes can be directly 
accessed using a Hanbury Brown- Twiss setup by looking 
at photon emission in a single output arm only, i.e. by 
measuring the g^ 2 ) of the signal or idler arm. In such an 
investigation however, jitter and arrival-time discretisa- 
tion smear the result. Due to the aforementioned dispro- 
portion between the time resolution of available single- 
photon detectors and the SPDC coherence time, most ex- 
perimental investigations have exploited fs-pulsed sources 
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and strong filtering to increase artificially the coherence 
time of the measured field. Positive evidence in single- 
photon regime was first reported by Tapster and Rarity 
who compared quantitative predictions based on a 
specific experimental hardware with measurements per- 
formed with various filters; with the narrowest one, a 
very convincing bunching peak of 1.85 could be demon- 
strated. Six years later, dc Ricdmatten et al. [3 (al- 
most) replicated the experiment, with similar shortcom- 
ings and results, concluding that the amount of bunching 
depends only on the ratio t c /t p . In both cases an appar- 
ent transition from thermal to Poissonian statistics was 
observed when r c was varied below r p . Soderholm et al. 
fl7l ] studied the dependence of single- and double-click 
rates on pump power in the two cases r c ~ 1.75r p and 
r c ~ 0.6r p , and found good agreement with a thermal 
model in the first case, and an intermediate behaviour 
in the second case. On the other hand, Mori et al. [l8[ 
studied a source with very long pump pulses (40ps) and 
broad filtering leading to a rather short coherence time r c 
(~ 120fs), and failed to detect any significant bunching. 
Similarly, in a recent article, Avenhaus et al. 15], using 
60ps laser pulses, observed Poissonian statistics for the 
SPDC process (r c not stated). 

In an experiment where, contrary to the single-photon 
regime, several photons per pulse are created, Paleari 
et al. [2CJ] were able to demonstrate, with careful fitting of 
the measured distribution, a thermal behaviour, although 
their ratio t p /t c was rather unfavourable (in the range 
10-20). Vasilyev et al. (2l[ used a strongly filtered opti- 
cal homodyne tomography to extract the photon-number 
distribution of a SPDC process and compare it with a 
thermal one, finding almost perfect agreement [32j . 

SPDC sources have also been implemented using 
continuous-wave (cw) pumping; this alternative approach 
is still relatively unexplored due to the technical diffi- 
culty of measuring very short correlation times. Larchuk 
et al. [IH tried to investigate sources in the single-photon 
regime with a single detector and delay lines but failed 
to obtain any evidence because the dead time, about 
lfis, makes the apparatus blind to the interesting re- 
gion t < lps. Super-Poissonian behaviour was found by 
Zhang et al. [22| by directly observing the photocurrents 
of orthogonally polarised twin beams from a continuously 
driven KTP crystal in a cavity. 

Measurements of the marginal, i.e. one-arm, SPDC 
field are not to be confused with measurements of the 
better known statistics of the whole bi-photon state, as 
well as conditioned statistics, which have been experi- 
mentally confirmed, among others, by [23l. |24l [25l. 126] . 



C. Quantum key distribution (QKD) and SPDC 
bunching 

Sources based on continuously pumped SPDC are of- 
ten employed in quantum cryptographic devices, gener- 
ating the multi-photon pair state of Eq[3J Correlations 



between multi-photon pairs constitute a potential threat 
to the security of QKD as the exchange is then suscep- 
tible to a photon- number-splitting attack [27]]. The case 
of entangled-photon sources [28| is still relatively unex- 
plored, and its solution depends on a thorough under- 
standing of the extent to which multiple SPDC pairs are 
correlated. 

In entangled-photon sources, the state, to first or- 
der, corresponds to a pair completely entangled in the 
information-carrying degree of freedom (polarisation); 
e.g. for typc-II SPDC |^~) = QH„Vl) - \V B ,Hi))/y/2. 
When the second order of the expansion (Eq(3|) is con- 
sidered (two pairs), two limiting cases emerge: sister 
pairs are in well distinguishable temporal modes, and are 
therefore completely uncorrelated in polarisation, with 
all polarisation combinations equally likely, 

[sisters) J^ + '^f -^.ra), (5) 

while twin pairs arc in the same temporal mode (within 
r c ), and 

. . \2H,2V) + \2V,2H)-\HV,VH) 
(twins) = -j= . (6) 

For twin states it has been shown that the probability 
that the two signal photons have the same polarisation 
is | @; therefore, whereas no kind of photon-number- 
splitting attack is possible with sister states, twin states 
are potentially dangerous for quantum key distribution. 
The situation is analogous for type-I SPDC. 

The works of Tsujino et al. [23] and Ou [3(| have shown 
that polarisation correlation and one-arm bunching arc 
two sides of the same phenomenon: the peculiar proper- 
ties of the state space of identical particles imply that, for 
delays shorter than r c , signal photons with the same po- 
larisation occur more frequently than expected for clas- 
sically uncorrelated objects (as differently polarised pho- 
tons). 

Summarising, the study of one-arm g^ 2 ) hi entangled- 
photon sources [33| gives information about the maxi- 
mum possible extent of the security threat posed to quan- 
tum key distribution by higher orders of SPDC. 



II. EXPERIMENTAL SETUP AND DATA 
ACQUISITION 

In our SPDC source [34| (FigJTJ) , a 532nm cw laser 
pumps a single 30mm temperature-stabilised nonlin- 
ear crystal (periodically-poled KTP) for type-I down- 
conversion (532nm — * 810nm + 1550nm). The signal 
and idler photons have fixed polarisations, central wave- 
lengths of respectively 810nm and 1550nm, and are emit- 
ted collinearly. The signal coherence time, as measured 
with a Michelson interferometer (Fig(3J), is r c ~ 2.8ps 
FWHM, corresponding to a bandwidth of less than lnm. 
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Signal and idler photons are separated according to 
their frequency, and coupled into single-mode optical fi- 
bres. Signal photons (at 810nm) are then recollimatcd 
and sent to a non-polarising balanced beam splitter. The 
resulting beams are collected into multi-mode fibres. The 
coupling efficiency is greater than 90% as the light is cou- 
pled from single-mode fibres into multi-mode ones. The 
multi-mode fibres direct the beams into a PerkinElmer 
SPCM-AQ4C single-photon detector array. The parame- 
ters of each detector are as follows: quantum efficiency 
at 810nm ~ 50%, dark-count rate < 500Hz, dead time 
~ 50ns, and detector saturation level above 1MHz. All 
detector clicks are recorded by a time-tagging unit (TTU) 
(TT8, smart systems Division, ARC) with a common 
time basis and an intrinsic resolution of t tt = 82.2ps. 
Time tags are then processed to define coincidences. 

One of the two fibres is a 100m spool, so that it oper- 
ates as a delay line of approximately 500ns (this delay is 
removed by software during data analysis). The purpose 
of the delay line is to make photons that happen to come 
with a picosecond delay to impinge on detectors at suffi- 
ciently different times to avoid electronic cross-talk and 
increase sensitivity. 

The timing jitter of the overall detection unit was char- 
acterised by measuring correlations in arrival times of 
signal-idler photon-pairs from an additional degenerate 
SPDC source at 810nm (405nm ->■ 810nm + 810nm). 
This auxiliary photon-pair source is based on a Sagnac 
configuration [3 11 ] - For the present measurement of jit- 
ter, it was operated to yield ~ 1MHz detected single 
rates, below detector saturation. The coincidence rate of 
the photon-pair source was approximately 100kHz. The 
combined jitter of the detection unit is estimated to be 
tj = 640ps ^> t c . This estimated FWHM is however not 
used in further experimental analysis, and no assumption 
on the shape of the jitter is needed. Instead, the fully 
normalised jitter curve, as seen in the cross-correlation 
in Fig|2J is used to characterise the detection unit, the 
only extracted parameter being the value of the area un- 
der the jitter curve (611ps). Several measurements of 
the jitter histogram were taken, and the deviation on the 
jitter area is estimated to be less than 3%. 

A cross-correlation histogram is obtained from record- 
ed detector clicks. A coincidence from detector D\ to 
detector D 2 with delay At = t 2 — t\ is counted if Di 
clicked at time t\, and D 2 clicked at time t 2 > t\, irre- 
spective of any other click. Coincidences with the role of 
detectors reversed are defined accordingly, and the two 
functions are joined by arbitrarily defining the delays At 
from D 2 to D\ as negative. 

The coincidence density for uncorrelatcd photons, e.g. 
when At is large with respect to any coherence time, 
is the product of the individual experimental detection 
rates Ai and A2 of D\ and D 2 (which include the effects 
of dead times and dark counts). Therefore, every bin of 
the coincidence histogram will contain, on average, and in 
the case of uncorrelatcd photons, C = \\\ 2 t tt T counts, 
where T is the duration of a data-taking run. Since 
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FIG. 1: (colour online) The experimental setup consists of a 
SPDC source of photon pairs, a balanced beam splitter and 
two single-photon silicon detectors, combined in a Hanbury 
Brown and Twiss 0] configuration. 



C contains both single rates Ai and A2, it is directly 
proportional to the square of the pump power. 



III. MODEL OF EXPECTED PHOTON 
BUNCHING AND MEASUREMENTS 

Since the jitter Tj of the detection unit is much larger 
than the coherence time r c of the down-converted light, 
any <7( 2 ) peak will be strongly smeared. A rough estimate 
for the residual peak height is given by the ratio of the 
coherence time to the jitter t c /t) ~4x 10 -3 . 

For a more accurate calculation, we model how a theo- 
retical thermal bunching peak {g/%) (0) = 2) will be trans- 
formed due to our instruments and then compare it with 
the experimental result. SPDC-bunching is characterised 
by the same relation between and the first order co- 
herence <7(i) that holds for chaotic light 0, namely: 



9(2){t) = l + |5(i)(r)| 2 



(7) 



The of the signal beam was measured in a Michel- 
son interferometer and can be seen in FigJ3] The setup for 
this measurement consisted of a Michelson interferometer 
with a motoriscd movable mirror and some limiting aper- 
tures in both arms to increase visibilities. The maximum 
visibility of the intcrferogram was measured to be 90%. 
The mirror is set to move at a constant speed of 0.002 
mm/s over a range of 3mm, and detection counts were 
recorded over integration windows of 2 ms. At unequal 
path lengths, the detector count rate at the output of the 
interferometer was in the range of 200kHz. The interfer- 
ence fringes in FigJ3]wcrc obtained from the raw data by 
subtracting the background rates (around 5kHz), and are 
displayed with an averaging of 30 data points to reduce 
fluctuations away from the centre of the interfcrogram. 

The envelope of the interferogram as shown in FigJ3J 
that is 13(1)1, is calculated from the visibility of the inter- 
ference fringes in the raw data (see inset in Fig(3|). An es- 
timate of the coherence time r c ~ 2.8ps is obtained from 
the graph. After rescaling the data points of the 13(1)1 
envelope to compensate for the reduced visibility (90%) 
in the interfcrogram, an area of 2.17ps was found for the 
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FIG. 2: (colour online) Second-order correlation function 
between the idler and signal arms of an auxiliary degener- 
ate SPDC source at 810nm (405nm 810nm + 810nm), 
rescaled to unity. Since the time correlation between pho- 
tons of the same pair is tight (sub-picoseconds), the FWHM, 
approximately 640ps, is essentially an estimate of the com- 
bined jitter of the two silicon detectors and the time-tagging 
unit. 



FIG. 3: (colour online) of the 810nm arm of the SPDC 
source as measured in a Michelson interferometer vs. optical 
path difference, gm was calculated from the fringe-visibility 
of the interferogram. The experimental visibility at zero opti- 
cal path difference is slightly less than 1. When the maximum 
of the visibility is rescaled to one, the total area under the 
squared envelope is ~ 2.2ps. The inset shows the visibility 
fringes around the central part of the interferogram. 



squared envelope |<7(i)| 2 - The error on this value can be 
estimated by considering the area of several interference 
scans in identical experimental conditions; a conserva- 
tive estimate is a 4% error on the calculated area for the 
squared envelope. 

To model the washed out <7( 2 ) peak, we note that the 
total number of coincidences in the excess peak is pre- 
served despite the jitter. Hence we expect to observe a 
smeared peak with the same temporal profile as that in 
FigEl but with an area equal to that of |ff(i)| 2 This 
method does not rely on any a priori knowledge of the 
actual shapes of <?(i) or g(2), but only on a ratio of ar- 
eas. The estimated coherence time r c and jitter Tj of 
the detection unit from the FWHM of their respective 
histogram are not used in the calculation, neither is any 
assumption on curves shapes necessary. 

The correlation function g^ 2 ) of the field from one arm 
of a SPDC source was measured using the Hnbury Brown- 
Twiss-like setup introduced in the previous section. A 
cross-correlation histogram for positive and negative de- 
lays was obtained from detector clicks and the bunching 
peak emerged in the region At ~ over the plateau of 
accidental coincidences. The average count number of 
the plateau was then used to normalise the histogram. 
Experimental data were collected during a 16 hours data- 
taking session with an average count rate of ~ 1 MHz in 
each detector (below saturation), giving a plateau count 
of N ~ 4.4 x 10 6 events per bin. The observed rates 
used for jitter characterisation were set to be identical 
to the rates in this measurement, since detector jitter 
is dependent on actual count rates. The double pair 
term in Eq[3] scales with the square of the pump power, 
as does the plateau (see the estimate C at the end of 



Sec|H]), hence the normalised g^2) function is indepen- 
dent of pump power. The same argument can be used to 
prove the peak to be independent of optical losses in the 
setup. 

In Fig. 31 we show a comparison of the measured data 
to the expected gr^ of our model. Data points corre- 
spond to the experimental g^) histogram of the marginal 
SPDC field. The solid line is the expectation for the 
bunching peak of a field with thermal statistics. This 
line was obtained by vertically rescaling the jitter func- 
tion shown in Figj2]to yield the same area as |ff(i)| 2 . The 
comparison does not involve any free parameter to be fit- 
ted: the theoretical model is only based on the specific 
and independently measured parameters of the setup. 

The peak was shifted laterally by ~ 50ps for better 
comparison with the experimental data; this is however 
less than the intrinsic resolution t tt . Statistical fluctua- 
tions of the number of events in a bin of the histogram 
are of the order of y/~N. The cumulative error of ±5% 
on the expected g^) curve, arising from a conservative 
error estimation from the jitter and |<?(i)| 2 data, is shown 
in FigH] as a shaded area. Even this 10% error margin 
fits the data points, and therefore the good match of the 
predicted peak is not accidental. 

The emerging bunching from the SPDC field is in excel- 
lent agreement with the theoretical estimation of a ther- 
mal statistics. The peak protrudes 6 standard deviations 
above the Poissonian limit of g^ 2 ) = 1 and is therefore 
experimentally confirmed. 
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FIG. 4: (colour online) Second-order correlation function 
(signal-signal <?(2)) of the 810nm arm of a SPDC source, nor- 
malised over a measured plateau of N counts in a bin. The 
solid line represents the washed-out bunching peak (original 
height of 2), obtained from the jitter plot rescaled to an area 
of 2.17ps. Error bars, representing statistical fluctuations of 
the measured data, are set to y/N. The shaded grey area, 
only visible at the very top of the peak, represents the uncer- 
tainty in the model prediction. The inset shows an enlarged 
section of the plateau. 

IV. CONCLUSIONS 

We experimentally measured the marginal temporal 
second-order correlation function 3(2) (t) of SPDC light 
with continuous wave excitation in the single-photon 
regime. To our knowledge, the results represent the first 
direct observation of photon bunching for such a field, 
and strongly confirm its thermal character. Even with an 
experimental excess peak g^ 2 ) (0) as low as 10 -3 , our ex- 
perimental accuracy is sufficient to properly observe the 
residual bunching peak, which can be modelled solely on 
the experimental timing jitter and the output of a first- 
order interferometric measurement (<7(i))- With much 



higher timing precision (e.g. lower detector jitter), it is 
expected that a full peak with the shape given by Eq. 
is recovered. The relevant parameter is the ratio between 
the coherence length t c and the measurement resolution 
(the jitter Tj in the cw case). 

We therefore conclude that thermal photon statistics 
of the marginal SPDC field can be observed both in the 
fs-pulsed and in the cw regime, the latter being more 
practical and more affordable. A multimode excitation 
does not change the statistics of SPDC but alters only 
the experimental conditions, so that without a careful 
analysis an apparent Poissonian character is observed (as 
presented in some works reviewed in section \TB\i , instead 
of the inherent thermal statistics. We have shown in this 
paper that even with far-from-ideal detection modules, 
and without any additional spectral filtering, we could 
demonstrate as a proof of principle that the true nature 
of the field remains accessible. 

It is hoped that the results presented here stimulate 
further investigations into the nature of SPDC light and 
its thermal properties. For cw pumping, there is no clear 
consensus on the actual state of multi-photon pairs, nor 
is its extension to entangled states produced by cw pump- 
ing well understood. Similar experiments would provide 
experimental data to confirm that bunching implies po- 
larisation correlations. In addition, cw pumping allows 
a continuous temporal investigation in a time domain 
not accessible to experiments with pulsed pumping and 
hence, with an improved timing resolution, access to the 
shape of g( 2 ) > an d not only to its integral (a single data 
point) as in pulsed setups. 
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